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■ The Rayleigh-Benard theory by Grossmann and Lohse [J. Fluid Mech. 

407, 27 (2000)] is extended towards very large Prandtl numbers Pr. The 
. Nusselt number Nu is found here to be independent of Pr. However, for fixed 

Rayleigh numbers Ra > f0 10 a maximum around Pr ~ 2 in the Nu(Pr)- 
dependence is predicted which is absent for lower Ra. We moreover offer 
the full functional dependences of Nu(Ra, Pr) and Re(Ra, Pr) within this 
extended theory, rather than only giving the limiting power laws as done in 
ref. Q. This enables us to more realistically describe the transitions between 
Q | the various scaling regimes, including their widths. 
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In thermal convection, the control parameters are the Rayleigh number Ra and the 
Prandtl number Pr. The system responds with the Nusselt number Nu (the dimensionless 
heat flux) and the Reynolds number Re (the dimensionless large scale velocity). The key 
question is to understand the dependences Nu(Ra,Pr) and Re(Ra, Pr). In experiments, 
traditionally the Prandtl number was more or less kept fixed [^-^]. However, the recent 
experiments in the vicinity of the critical point of helium gas |5|,|6[ and of SF$ or with 
various alcohols || allow to vary both Ra and Pr and thus to explore a larger domain of the 
Ra—Pr parameter space of Rayleigh-Benard (RB) convection, in particular that for Pr 3> 1. 
While the experiments of Steinberg's group |7j suggest a decreasing Nusselt number with 
increasing Pr, namely Nu = 0.22Ra°- 3±om Pr'°- 2±om in 10 9 < Ra < 10 14 and 1 < Pr < 93, 
the experiments of the Ahlers group suggest a saturation of Nu with increasing Pr for fixed 
Ra, at least up to Ra = 10 10 ||. The same saturation (at fixed Ra = 6 • 10 5 ) is found in the 
numerical simulations by Verzicco and Camussi [10] and Herring and Kerr [[II] . 
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The large Pr regime of the latest experiments has not been covered by the recent theory 
on thermal convection by Grossmann and Lohse (GL, 0), which otherwise does pretty well 
in accounting for various measurements. In particular, it explains the low Pr measurements 
of Cioni et al. [|J (Pr = 0.025), the low Pr numerics which reveal Nu ~ p r 014 for fixed Ra 
10|) TTf] , and the above mentioned experiments by Niemela et al. || and Xu et al. ||. 



In the present paper we extend the GL theory in a natural way to the regime of very 
large Pr, on which no statement has been made in the original paper Q. We find Nu 
to be independent of Pr in that regime. We in addition present the complete functional 
dependences Nu(Ra, Pr) and Re(Ra, Pr) within the GL theory, rather than only giving the 
limiting power laws and superpositions of those as was done in This enables us to more 
realistically describe the transitions between the various scaling regimes found already in 

!■ 

Approach: To make this paper selfcontained we very briefly recapitulate the key idea of 
the GL theory, which is to decompose in the volume averages of the energy dissipation rate 
e u and the thermal dissipation rate ee into their boundary layer (BL) and bulk contributions, 
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£u — £u,BL + €u,bulki (1) 
,bulk ■ (2) 

For the left hand sides the exact relations e u = j^{Nu — l)RaPr~ 2 and e e = KjjNu are 
used, where v is the kinematic viscosity, k the thermal diffusivity, L the height of the cell, 
and A the temperature difference between the bottom and the top plates. Next, the local 
dissipation rates in the BL and in the bulk (right hand sides of eqs. (P and (0)) are modelled 
as the corresponding energy input rates, i.e., in terms of U, A, and the widths X u and Xg 
of the kinetic and thermal boundary layers, respectively. For the thickness of the thermal 
BL we assume Xg = L/(2Nu) and for that of the kinetic one X u = Lj (Ay/Re) as it holds 
in Blasius type layers ]nj; as for the prefactor 1/4 cf. Sec. 4.3. For very large Ra the 
laminar BL will become turbulent and X u will show a stronger Re dependence. Note that 
whereas the thermal BLs only build up at the top and bottom wall, the kinetic BL occurs 
at all walls of the cell and therefore the contribution of € u bl to e u is larger than a simple 
minded argument would suggest. The two eqs. ([]]) and (0) then allow to calculate the two 
dependent variables Nu and Re as functions of the two independent ones Ra and Pr. 

Input rate modeling: The modeling of the dissipation rates on the rhs of eqs. (|l|) and (||) is 
guided by the Boussinesq equations. Depending on whether the BL or the bulk contributions 
are dominant, one gets different expressions on the rhs of eqs. ([!]), (|2|) and thus different 
relations for Nu, Re vs. Ra, Pr, defining different main regimes, see ref. [0] and fig. |I[ 

The rhs thermal dissipation rates depend on whether the kinetic BL of thickness X u is 
within the thermal BL of thickness Xg (X u < Xg, small Pr) or vice versa (X u > Xg, large 
Pr). The line X u = Xg, corresponding to Nu = 1\fRe~, splits the phase diagram into a lower 
(small Pr) and an upper (large Pr) part, which we shall label by "f and 'V. 

We first consider X u < Xg (i.e., small Pr, regime "f ). Then (see M) 
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eu, BL - ~ (4) 

UA 2 A 2 

ee,buik 7 n—PrRe, (5) 

ee,BL~ K^(RePr) 1 / 2 . (6) 

The last expression is concluded [jl), j3| [l~3|| from the heat transfer equation u x d x 9 + u z d z 9 = 
nd 2 9, which implies U/L ~ k/Xq giving Re^^Pr 1 / 2 ~ Nu. 

If now larger Pr are considered, the kinetic boundary will eventually exceed the thermal 
one, X u > Xg, upper range "m". The relevant velocity at the edge between the thermal BL 
and the thermal bulk now is less than U , about UXg/X u . To describe the transition from X u 
being smaller to being larger than Xg we introduce the function f(x) = (1 + x n )~ x l n of the 
variable xg = X u /Xg = cgNu/2\^Re, f being 1 in the lower range "f (small Pr) and 1/xg 
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FIG. 1. Phase diagram in the Ra — Pr plane, with data points included where Nu has been 
measured or numerically calculated. Squares show measuring points of Chavanne et al. ||, dia- 
monds those by Cioni et al. Q, circles those by Niemela et al. ||, the very thick lines those by 
Xu et al. ||, and the triangles are those points for which Verzicco and Camussi did full numerical 
simulations [1C|. The long-dashed line is the line X u = Xg. The thin dotted line denotes where the 
laminar kinetic BL becomes turbulent. As extensively discussed in ref. [JD], the exact onset of this 
instability strongly depends on the prefactors used when calculating this type of phase diagram. 
For this phase diagram the same prefactors as in ref. O have been chosen. 



in "w" (large Pr), respectively, eg is about 1. The relevant velocity then is Uf(xg). We take 
n = 4 to characterize the sharpness of the transition. 
This generalizes (5), (6) to 

A 2 

€g M i k ~ K—PrRef(c e Nu/2VRe), (7) 



tg, BL ~ k— yf PrRef (cgNu/2^) , (8) 

while e U)bu i k and e U)B L are still given by (3), (4). Introducing (7), (8), (3), (4) into (1),(2) 
leads to the Ra,Pr dependences of Nu,Re in the upper regime "it". The pure power laws 
Nu(Ra, Pr) and Re(Ra, Pr) in both the lower (X u < Xg, small Pr) and the upper "w" 
(X u > Xg, large Pr) regimes are summarized in table I. 

Very large Pr regime: We now extend the theory to very large Prandtl numbers. As 
long as Ra > Ra c = 1708 there still is wind, even for very large Pr, as Ra c is independent 
of Pr. However, to keep Ra fixed, one has to increase the temperature difference A with 
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TABLE I. The power laws for Nu and Re of the presented theory. The regimes II U and 
are not included as they do not or hardly exist for the choice of prefactors. 

increasing v or Pr. Pr ^> 1 will lead to a smaller and smaller large scale wind Re. The flow 
will eventually become laminar throughout the cell. X u can no longer continue to increase 
according to \ u ~ Re~ l l 2 with decreasing Re, but will saturate to a constant value of order 
L. This is the central new point of the very large Pr regime. In |L| we had assumed that 
this happens at about Re = 50. This corresponds to a saturation at \ u = L/A^/Re = L/28. 
To model the smooth transition to the very large Pr regime beyond the line Re = 50, i.e., 
if X u = Lf Ay/Re approaches X u = L/28 we use the crossover function g{x) = x(l + x n )^ 1 ^ n 
of the crossover variable xl = 28X U / L = 7/ \f~Re~, and again n = 4. The function g increases 
linearly, g{x£) = xl, below the transition (x^ small) and is 1 in the very large Pr regime 
with Re < 50. In the above modelings for the local dissipations rates we have to replace 
each A„ by g(x L )L/28. 

The resulting formulae, given momentarily, will lead, depending on Ra, to three new 
regimes, valid for very large Pr, denoted as 1^, and 1 11^, see Fig. [I] and Table I. 

While eq. (3) for e u ^ u ik is assumed to still hold in the very large Pr range (where e u ,buik 
of course hardly contributes to e u due to the large extension of the kinetic BLs), eqs. (4), 
(5), and (6) have to be generalized. First generalize (4) for € U) bl, 

U 2 v 3 Re 2 

Cu,BL ~ V— — — ~ — — (9) 

g{x L )L 2 L 4 g(7/VRe) 

Next ee,buik- Above the wind velocity U in (5), which sets the time scale of the stirring, has 
already been generalized to Uf(X u /\e). This equals U itself in "f and UXe/X u in the "it" 
regimes. Now in addition the explicit A„ is to be replaced by g(xL)L/28, i.e., 

A 2 / Nu ( 7 \\ 

ee Ml k~K-PrRe f\—g\-—^. (10) 

This simplifies for large enough Ra (therefore large /-argument) and very large Pr (thus 
large g-argument) to 

A 2 PrRe /11X 

ee,bui k ~Kj2 -jy— • (11) 
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Inserting (9) and (11) into the rhs of (1) and (2) leads to the new power laws describing the 
heat flux and the wind velocity in the regime II loo beyond III U , cf. Figure 1, 



Nu ~ Ra 1/3 Pr°, Re ~ Ra 2/3 Pr 



(12) 



Finally ee,BL'- in the thermal boundary layer range beyond Ig, relevant for medium Ra, 
eq.(6) stays valid, because its derivation did not involve X u and also / = 1. The range 
has to be described by eqs.(l),(2),(9) (with g = 1), and (6), resulting in the same power 
laws as in regime IIIoo, be., eqs. (|T2]). For Pr values above regime I u eq.(6) no longer holds. 
It originated from the heat transport equation. There we have to use now Uf(xg) instead 
of merely U. The balance from the heat transfer equation then reads Uf(xg)/L ~ k/A 2 ,. In 
the /-argument xq = X u / Xg the kinetic BL width X u and therefore the crossover function g 
appears, leading to 



PrRef ( — g[ 

1 U y WRe' 




Nu 2 . 



(13) 



In the very large Pr regime (where g(xi) = 1) above I u (where f(xg) = Xg 1 ) one obtains 
from flT3"| ) the relation (RePr) 1 / 3 ~ Nu, valid in I>. Together with (1), (2), and (9) one 
derives the new scaling laws in the interior of II 
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Nu ~ Ra 1/5 Pr°, Re ~ Ra 3/b Pr 



(14) 



The scaling behavior Nu ~ Ra 1 / 5 has earlier been suggested by Roberts [Q. Note that in 
all three very large Pr regimes Nu does not depend on Pr. Furthermore the Constantin- 
Doering [13] upper bound Nu < constRa 1 ^ 3 ^! + log Ra) 2 / 3 , holding in the limit Pr — > oo, 
is strictly fulfilled. 
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FIG. 2. (a) Nu as a function of Pr according to theory for Ra = 10 6 , Ra = 10 7 , ...to 
Ra = 10 16 , bottom to top. (b) Nu/Ra 1 ^ as function of Ra for Pr = 10" 4 , Pr = 1(T 3 , Pr = 10~ 2 , 
Pr = 10 _1 , and Pr = 10° (solid lines, bottom to top) and for Pr = 10 1 , Pr = 10 2 , Pr = 10 3 , and 
Pr = 10 4 (dashed lines, top to bottom). 
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Nu and Re in the whole parameter plane: Plugging now the general expressions for the 
local dissipation rates (3), (9), (10), and ([TJ) into the balance eqs. ([!]) and (|2]) finally results 
in 

Re 2 

NuRaPr~ 2 = Cl =- + c 2 Re 3 , (15) 

g(7/VRe) 



Nu = c 3 Re 1/2 Pr 1/2 

Here we have added dimensionless prefactors where appropriate to complete the modeling 
of the dissipation rates. In ref. M c\ through c 5 were already adopted to Chavanne et 
al.'s experimental data ||. The result was c\ = 1028, c<i = 9.38, c 3 = 1.42, C4 = 0.0123, 
1/cg = C5 = 2.0, see eqs. (4.4) through (4.9) of Note that the q may depend on the 
aspect ratio and are not universal. 

The set of eqs.(|i~5"D and fll6|) is the second main result of this paper. It allows to calculate 
Nu(Ra, Pr) and Re(Ra, Pr) in the whole Ra — Pr parameter space, including all crossovers 
from any regime to any neighboring one. Technically, eq. ( |l"5f) is solved for Nu(Ra, Re, Pr) 
which then is inserted into eq. (|T6|), leading to one implicit equation for Re(Ra, Pr). 

All limiting, pure scaling regimes which can be derived from eqs. ( |i~5|) and (|16|) are listed 
in table I which extends table II of ref. fflj, now including the three new regimes J^, 1^, and 
1 1 loo. The corresponding phase diagram is shown in figure [I], completing that of |l[ towards 
very large Pr. 

Though in the phase diagram we have drawn lines to indicate transitions between the 
regimes, defined by either € u ,bl = e u ,buik or e e< BL = £e,buik or ^« = ^0 or ^« = L /28, we note 
that the crossovers are nothing at all but sharp. All transitions are smeared out over broad 
ranges, the more, the more similar the scaling exponents of the neighboring regimes are. 

Discussion of Nu(Pr) and Nu(Ra): The functions Nu(Pr) (for fixed values of Ra) and 
Nu(Ra)/ Ra l / A (for fixed values of Pr) resulting from eqs. (|T5"|) - (|T6| ) are shown in figure |2|. 
Indeed, for Ra ~ 10 6 — 10 9 the Nusselt number saturates with increasing Pr and the regime 
I u with Nu ~ p r -V 12 j s suppressed, in full agreement with the experimental and numerical 
results P|-4TT|. However, for larger Ra beyond ~ 10 11 the theory predicts a maximum for 
the curve Nu(Pr). The decreasing branches of the curve for Ra ~ 10 12 and Pr w 50 may 
be consistent with above mentioned experimental results by Ashkenazi and Steinberg [0. 
This observation may resolve the apparent descrepancy between the Ahlers et al. and the 
Steinberg et al. data: Simply different regimes in the Ra — Pr phase space are probed. 

The importance of transition ranges is highlighted in fig. 0b. E.g., in regime with 
Nu ~ Ra 1 / 5 (for fixed Pr) the corresponding scaling exponent 1/5 — 1/4 = —1/20 only 
becomes observable for very small Pr w 10~ 3 — 10~ 4 . At Pr = 10~ 2 the roughly four 
decades of regime /// (see figure 1) between regime // and IVi, which both have a stronger 
Ra dependence of Nu, are not sufficient to reveal the scaling exponent. And at Pr = 10 _1 
the roughly 2.5 decades of regime 11\ are nowhere sufficient to lead to a local scaling exponent 
dNu/dRa smaller than 1/4! 

Similarly, for Pr = 1 only for very large Ra > 10 15 pure scaling Nu ~ Ra 1 / 3 is revealed. 
For smaller Ra regimes I u and I\ with Nu ~ Ra 1 / 4 strongly contribute, resulting in an 
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FIG. 3. Nu(Ra, Pr(Ra)) / 'Ra 1 / 4 along the curve Pr(Ra) given by the experimental restrictions 
of either the Niemela et al. experiment || (boxes in the inset and solid line in the main figure; the 
experimental data itself are shown as boxes) or the Chavanne et al. experiment |5| (circles in the 
inset and dashed line in the main figure; the experimental data itself are shown as circles). The 
simple parameterizations on which the curves in the main figures are based are also indicated. 

effective local scaling exponent (increasing with Ra) in the range between 0.28 ("2/7") and 
0.31, just as observed in experiment f|,|3|,|^, |]. 

Direct comparison to experiment: Finally, let us point out that the present approach 
simplifies the comparison with experimental data: In experiment, it is hard to vary either 
Ra or Pr over many decades and at the same time to keep the other variable fixed. So 
most measurements are done along curved lines in the phase space figure [I], mixing the Ra 
and Pr dependences. Now eqs. ([15]) and ([U|) allow to calculate Nu and Re along such 
a curve Pr(Ra) given by the experimental restrictions, which connect Ra and Pr. E.g., 
figure |3| shows Nu(Ra, Pr(Ra)) with Pr(Ra) as resulting from the Niemela et al. || and 
the Chavanne et al. || experiments. 

The present theory suggests that beyond Ra ~ 10 12 the Nusselt number in the Chavanne 
et al. experiment || should be smaller than in the Niemela et al. experiment || , in contrast 
to what is found. This result is also in contrast to our former speculation of ref. that the 
strong Nu ~ Ra 3 / 7 dependence of regime III M , to which the Chavanne et al. data are closer, 
would account for their experimental findings. What we had overseen is that the larger Pr 
numbers of the Chavanne et al. experiment (for fixed Ra) and the Nu ~ Pr -1 / 7 dependence 
of regime III U overcompensates the strong Ra dependence of Nu for the Chavanne et al. 
data close to regime II I u . 

This present theory cannot resolve the paradox between the Chavanne et al. and the 
Niemela et al. data. Possibly, different temperature boundary conditions have been applied. 
Possibly, Chavanne et al. have already observed the transition from a laminar to a turbulent 
kinetic BL (dotted line in the phase diagram fig. [I]). 
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